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Photoassociation with short laser pulses has been proposed as a technique to create ultracold 
ground state molecules. A broad-band excitation seems the natural choice to drive the series of 
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excitation and deexcitation steps required to form a molecule in its vibronic ground state from 
two scattering atoms. First attempts at femtosecond photoassociation were, however, hampered 
by the requirement to eliminate the atomic excitation leading to trap depletion. On the other 
hand, molecular levels very close to the atomic transition are to be excited. The broad bandwidth 



of a femtosecond laser then appears to be rather an obstacle. To overcome the ostensible conflict 
of driving a narrow transition by a broad-band laser, we suggest a two-photon photoassociation 
scheme. In the weak-field regime, a spectral phase pattern can be employed to eliminate the 
atomic line. When the excitation is carried out by more than one photon, different pathways in 
the field can be interfered constructively or destructively. In the strong-field regime, a temporal 
£>y phase can be applied to control dynamic Stark shifts. The atomic transition is suppressed by 

choosing a phase which keeps the levels out of resonance. We derive analytical solutions for atomic 



two-photon dark states in both the weak-field and strong-field regime. Two-photon excitation may 
thus pave the way toward coherent control of photoassociation. Ultimately, the success of such a 
scheme will depend on the details of the excited electronic states and transition dipole moments. 
We explore the possibility of two-photon femtosecond photoassociation for alkali and alkaline-earth 
metal dimers and present a detailed study for the example of calcium. 



I. Introduction 

Coherent control was conceived to solve the problem of optimizing the outcome of a chemical 
reaction^. Initial work was devoted to photodissociation^. In this study, we address the inverse 
problem of controlling the free-to-bound transition in photoassociation^ where atoms are assembled 
to a molecule by laser light. The underlying principle of coherent control is interference of different 
pathways that lead from the initial state to the final outcome. Typically, in a free-to-bound 
transition, the relative phase between the reactants is not defined^: If the wavefunction of the 



initial state can be written as a product of the wavefunctions of the atoms, the outcome of the 
binary reaction cannot be controlled. At very low temperatures, these considerations do not hold. 
Ultracold collisions are characterized by threshold effects^. The kinetic energy is very small, and 
the inner part of the scattering wavefunction is dominated by the potential energy. As a result, the 
colliding pair becomes entangled. Ultracold atoms are therefore the best candidates for coherent 
control of a binary reaction. The formation of molecules in ultracold atomic gases is presently 
attracting significant interest^. 

Previous proposals for short-pulse photoassociation were all based on one-photon transitions. 
Initially the goal was to optimize the excitation process in photoassociation by introducing a 
chirp 10,11 . In order to produce molecules in their ground electronic state, a pump-dump scheme 
was studie d 12 ' 13114 : A first pulse excites two atoms and creates a molecular wave packet. This is 
followed by free dynamics on the electronically excited state bringing the wave packet to shorter 
internuclear distances. There a second pulse stabilizes the product by transferring amplitude to the 
ground electronic state. The pump-dump scheme can be supplemented by a third field to engineer 
the excited state wave packet dynamics to create favorable conditions for the dump pulse 15 . 

Experiments aiming at femtosecond photoassociation with a one-photon near-resonant excita- 
tion were faced with the obstacle that weakly bound molecular levels are very close to an atomic 
resonance. The challenge is to populate these levels without exciting the atomic transition. The 
weakly bound molecular levels possess the biggest free-bound Franck-Condon factors. However, 
the atomic transition matrix elements are several orders of magnitude larger. Excitation of the 
atomic transition is followed by spontaneous emission and trap loss, i.e. it leads to a depletion of 
the trap 16 ' 17 . Pulses with bandwidths of a few wavenumbers corresponding to transform-limited 
pulse durations of a few picoseconds were proposed as a remed y 12 ' 13 . However, for such narrow 
bandwidths, pulse shaping capabilities have yet to be developed. Experimenters employing fem- 
tosecond pulses with bandwidths of the order of 100 cm -1 have resorted to suppressing the spectral 
amplitude at the atomic resonance frequency by placing a razor knife into the Fourier plane of the 
pulse shape r 16 ! 17 ' 18 . Due to the finite spectral resolution, this leads also to the suppression of spec- 
tral amplitude which could excite the molecular levels with the largest free-bound Franck-Condon 
factors as illustrated in Fig. [T] (right-hand side). On the other hand, the spectral amplitude from 
the peak of the spectrum as well as from just below the cut finds almost no ground state popula- 
tion to excite since the probability of finding two colliding atoms at short internuclear distance is 
extremely small (left-hand side of Fig. [T]). As a result, only a tiny part of the spectral amplitude 
contributes to the photoassociation signal 18 . 



Here we suggest to employ two-photon transitions for femtosecond photoassociation. The main 
idea consists in rendering the atomic transition dark by applying spectral or temporal phase control 
while exciting population into weakly bound molecular levels by two-photon transitions. The spec- 
tral resolution is determined by that of the pulse shaper, typically of the order of one wavenumber. 
This should allow for the population of weakly bound molecular levels with the best free-bound 
Franck-Condon factors. 

Two-photon control is based on constructive or destructive interference of all two-photon path- 
ways adding up to the two-photon transition frequenc y 19 ' 20 . It allows for the excitation of a very 
narrow atomic transition by a broad band femtosecond pulse^ 9 -. In particular, an anti-symmetric 
spectral phase optimizes non-resonant two-photon population transfer while a symmetric spectral 
phase yields a dark pulsed. For weak fields, the two-photon absorption yield can be calculated 
within perturbation theory, allowing for an analytical derivation of suitable phase functions. At 
intermediate fields, four-photon pathways additionally contribute to the two-photon absorption, re- 
quiring a higher order in the perturbation treatmen t 21 ' 22 . In both the weak-field and intermediate- 
field regime, rational pulse shaping is based on frequency domain argument a 19 ' 20121 ' 22 . For strong 
fields, a time-domain picture becomes more adequate since the dynamics are dominated by time- 
dependent Stark shift o 23i24i2S . The atom accumulates a phase due to the dynamic Stark shift. A 
suitable control strategy consists in eliminating this phase by applying a time-dependent phase 
function, i.e. by chirping the pulse 2 ^. This allows for maintaining a 7r-pulse condition despite the 
levels being strongly Stark shifted. 

These results for coherent control of two-photon absorption in atomic systems serve as our 
starting point to derive pulses which are dark at the atomic transition but optimize two-photon 
photoassociation, i.e. population transfer into molecular levels. The shape of the pulses yielding 
a dark resonance of the atoms are derived analytically: In the weak-field and intermediate-field 
regimes, the two-photon spectrum is required to be zero at the atomic transition frequency, while 
for strong fields a time-dependent phase allows to maintain an effective 2-7r-pulse condition for the 
atoms. The shaped pulses are tested numerically for molecule formation. The paper is organized as 
follows: An effective two-state Hamiltonian is derived in Section [U by adiabatically eliminating all 
off-resonant levels and invoking the two-photon rotating- wave approximation (RWA). In Section llTH 
conditions for pulses that are dark at the two-photon atomic transition are derived analytically. 
Photoassociation is studied in Sections HV] and |V] with Section |IV] reviewing possible two-photon 
excitation schemes for alkali and alkaline-earth dimers and Section [V] presenting the numerical 
study of two-photon photoassociation for calcium. Section |VT] concludes. 



II. Theoretical framework for multi-photon excitations 



We extend the treatment of multi-photon transitions in atom a 19i20i2:Li22123125 to include the vibra- 
tional degree of freedom, R, of a diatomic molecule. To this end, we write the nuclear Hamiltonian 
for each electronic state, = T + Vi(R) + u>i (setting h = 1), i.e. all potentials, VJ(R), go to 
zero asymptotically and the atomic excitation energies are contained in Wj. T denotes the vibra- 
tional kinetic energy. Assuming the dipole approximation for the matter-field interaction, the total 
Hamiltonian is given by 

H = H + p,E(t), (1) 

where fi denotes the dipole operator and E{t) the laser field with envelope Sit), E(t) = 
^EoS(t) (e2 v ^ e~ luJLt + c.c^j. For weak and intermediate fields, the dynamics under the Hamil- 
tonian, H, can be solved by perturbation theory. To lowest order, the two-photon absorption is 
given in terms of the two-photon spectral amplitude on resonanc o 19l2 ° l21 , 



So 
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E{(jj ge /2 + uj)E{u ge /2-oj)duj , (2) 



where ui ge is the two-photon resonance frequency and E(lo) the Fourier transform of E(t). For an 
atomic system, conditions for a dark two-photon resonance are derived analytically by setting the 
integral to zero, see below in Section IIII Al 

For strong fields, we follow the derivation of Trallero-Herrero et. a l 23 ' 25 . The matter Hamiltonian 
is expanded in the electronic basis, 



H = J2(f + Vi(k)+Uij 



(3) 



and the electronic states are separated into initial ground state, final excited state and intermediate 
states. If the one-photon detunings of the intermediate states, ^ m ,g/e ~ ^L, are large compared 
to the spectral bandwidth of the laser pulse, Au>l, the intermediate states can be adiabatically 
eliminated^ 3 -. Since the electronic excitation energies are much larger than the vibrational energies, 
the adiabatic elimination in the molecular model proceeds equivalently to the atomic case. Invoking 
a two-photon RWA, an effective two-state Hamiltonian is obtained, 

A(t)= ft + V g (R) + ^(t) X(t)e~^ \ 

\ x(t)e^) f + y e (R) + A 2P + cuf(t)/ ' 

where u>^ e (t) denotes the dynamic Stark shift of the ground and excited state, x(t) the two-photon 
coupling and A 2 p the two-photon detuning, A 2 p = uj eg — 2lol. The two-photon coupling is given 



in terms of the one-photon couplings and one-photon detunings of the ground and excited state to 
the intermediate levels, 

X {t) = -\El\S(t)\^ J ^^- (5) 

m m 9 L 

Similarly, the dynamic Stark shifts are obtained as 

^) = 4^I^)| 2 EW 2 ^S» i = 9,e. (6) 

m mi L 

The transition dipole matrix elements, fcj, and the transitions frequencies, uj m i, are in general in- 
dependent. They tend to a constant value (the atomic dipole moments and transition frequencies) 
at large internuclear distances that are important for photoassociation. The -R-dependence is 
therefore neglected in the following. 

Eqs. ([5]) and © contain couplings between the intermediate and the ground and excited states, 
but couplings between intermediate states are neglected. For very strong fields, these higher 
order couplings should be included. This type of adiabatic elimination is valid if the Stark shift 
is small relative to the energy spacing between the ground and excited states. The validity of 
the conditions for adiabatic elimination and two-photon RWA are easily checked in the atomic 
case by comparing the dynamics under the effective Hamiltonian, Eq. (jl]), to those under the full 
Hamiltonian, Eq. ([I]), including all intermediate levels. Possibly the Born-Oppenheimer potential 
energy surfaces in Eq. ([3]) show crossings or avoided crossings in the range of R that is addressed 
by the pulse. In such a case, Eq. (|2|) needs to be extended to include all coupled electronic states. 

For strong-field control, it is instructive to choose the rotating frame such that the phase terms 
of the coupling in Eq. (jl]) are expressed as time-dependent frequencies * and to time-dependently 
shift the origin of energy 23 . The transformed Hamiltonian is written, 

(t + Vq(R) x(t) \ 

H(t) = { 9 , (7) 

\ x (t) t + V e (R) + {S^(t) + A 2P + <p)J 

where the differential Stark shift 

5 s Ut)=u J s e (t)-oj s g (t) (8) 

is introduced. In the atomic case, strong-field control consists in locking the phase corresponding 
to the term in parenthesis in Eq. ([7]), i.e. in keeping 

*(*)= f €(r)dr + A 2P t + <p(t) (9) 

J — oo 



This corresponds to the fact that in the RWA a chirp can be expressed either as a phase term of the electric field 
envelope or as instantaneous laser frequency. 



constant. This can be achieved by proper choice of phase of the laser pulse, f(t), and it corresponds 
to maintaining resonance by instantaneously correcting for the dynamic Stark shifts. 



III. Two-photon dark pulses 

A. Solution in the weak- and intermediate-field regime 

At low field intensities, the population transfer is dominated by the resonant two-photon path- 
ways. Therefore control strategies are best understood in the frequency domain 21 . Starting from 
the Hamiltonian in the two-photon RWA, Eq. @, the condition for a dark pulse at the atomic 
transition can be derived analytically. The atomic transition corresponds to the asymptotic limit 
of Eq. ([3]), i.e. to a two-level system. Since the two-photon coupling, is proportional to the 
field amplitude squared, EQ\S(t)\ 2 , we look for a spectral phase such that the Fourier transform of 

/+oo r+oo 
S 2 (t)e lujt dt = / S(u/)S(w - J)dJ , 
-oo J — oo 

becomes zero at the two-photon resonance, F{u eg ) = (S(co) is the Fourier transform of the 
envelope function S(t)). We assume a two-photon resonant pulse, ujl = w eg /2. Then 

/+oo 
S(uj')S(2uj l - lo')cLu' (10) 
-oo 

/+oo 
S{u)l + uj)S(u!l — uS)duj . 
-oo 

The condition in Eq. (|10p is formally equivalent to that obtained in second order perturbation 
theor y 19 ' 21 . However, here we have derived it by considering the effective two-photon Hamiltonian, 
Eq. (jl]). The condition of Eq. ()10p is sufficient for a two-photon dark pulse as long as the states 
are not strongly Stark-shifted, i.e. it holds also at intermediate field intensities. This finding is 
rationalized by considering the next order contribution. The fourth-order terms are classified into 
on-resonant and near-resonant terms 2 ^. The dominant fourth-order term is the one that includes 
the resonant second-order term and is therefore canceled by the same condition. 

We assume the pulse shape to be Gaussian. Then S{uj) = A$e ^2 e l< t>{^) with the spectral 
bandwidth, and the two-photon Fourier transform becomes 

F(2co L ) = Al e^e^^+^-^W 

J — oo 

We now need to choose the spectral phase 4>(cu) such as to render F(2lul) zero. This can be achieved 
by a step of the phase 4>(uj) by 7 r 19 i 20 . 



We derive here the position of this step, u s . Let us assume that uj s > ojl, and denote Auj s 
uj s — u>l- Then 
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uj > Aw s , 




UJ L + Uj) + <p(0JL — U>) = < 
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uj < —Auj s , 


(11) 
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if 


— Auj s < uj < Auj s . 





The 7r-step implies a change of sign in the integral, and the dark pulse condition is rewritten 
accordingly, 

e ^duj+ / e ^duj- / e ^dw, 

-oo J Au s J—Auj a 

-oo _ ul 2 rAuj s _ w 2 

We recognize the error function up to a factor and obtain 



which determines Au; s and hence the step position, uj s . That is, if one chooses Auj s such that 
Erf(crAu s ) = 1/2, the two-photon spectrum is zero at the two-photon transition frequency, 
F(2ujl) = 0, and a dark pulse is obtained. Interestingly we can also combine the conditions 
of a symmetrical function and a phase step. This is achieved by defining two phase steps sym- 
metrical around ujl, i w s . Following a derivation along the same lines as outlined above, the 
identical condition for u;, is obtained. 



B. Solution in the strong-field regime 



At high intensities, the effective Hamiltonian varies during the pulse. A time-domain picture is 
therefore best adapted to develop control strategies. Population inversion in the effective Hamil- 
tonian can be achieved by adjusting the temporal phase of the laser such as to maintain a -zr-pulse 
conditio n 23 ' 24 . Here we turn this strategy upside-down to achieve a dark pulse, i.e. we ask for the 
27r-pulse condition. We restrict the notation to the atomic case for simplicity. The Hamiltonian, 
Eq. ([7]), can be transformed such that the diagonal elements are zerc- 2 ^ 

x(*)e~ i0(t) ' 
v x(*)e i$(i) 

Evaluating Rabi's formula for the final state amplitude leads to the condition 

*/ 



H(t) 



(13) 



xity^dt 



2kir. 



0,1,... 



(14) 



for zero population transfer to the excited state. We again assume a two-photon resonant pulse, 
uj ge = 2ujl. The condition, Eq. ([HP , is fulfilled given the phase <fr(i) obeys 

if t - At x < t < t + At x , 
$(t) = { (15) 
elsewhere . 

This corresponds to the following temporal laser phase 

.,. , -^-Io S 5(r)dT if t - At x < t < to + At x 
<PW = \ , (16) 
J ^(r)dr elsewhere 

For a Gaussian pulse, the time-dependence of the differential Stark shift is also Gaussian, and (pit) 
can be evaluated analytically. Inserting the definition of 3>(i), Eq. Q, into the condition, Eq. (|14|) . 
and evaluating the integrals yields a condition for At x , 



A pulse with the temporal phase defined by Eqs. (|16l) and (fTTj) suppresses population transfer by 
keeping the two levels out of resonance. 

Such a solution can also be obtained by applying local control^ to the effective Hamiltonian, 
Eq. (I13p . Local control theory allows to carry the analysis over to the full molecular problem. We 
can seek conditions which eliminate atomic two-photon transitions while maximizing the molecular 
population transfer. To this end, we define an operator to be minimized, 



at 



oo 

dR\e){e\, (18) 
Ro 



which projects onto the atomic levels, |e) denotes the excited electronic state. Ro is a distance which 
is larger than the bond length of the last bound level of the excited state potential. The objective 
to be maximized is defined by the projection onto the molecular levels for photoassociation, 



mol 



where \ipf) denote vibrational eigenstates of the electronically excited state. In local control theory, 
we seek the phase of the field such that 

J t (Pat) = j t (Pmol}>0, (20) 

i.e. atomic transitions are suppressed while the molecular population is monotonically increasing-- . 
Eq. (fTBT) is a specific analytical example of local control. The first part of the pulse up to time 



t < tQ — At x builds up excited state population and a phase relation between the ground and excited 
levels. At this instant in time the phase of the field jumps by it so that it becomes perpendicular to 
the phase of the instantaneous transition dipole, thus eliminating additional population transfer. 
The last part of the pulse restores the original population. 

IV. Two-photon excitation schemes for photoassociation 

Our goal is to populate molecular levels in a long-range potential while suppressing excitation 
of the atomic transition. We will restrict the discussion to homonuclear dimers where potentials 
with 1/B? behavior occur below the S + D asymptote* . In heteronuclear dimers, the potentials 
go as 1/i? 6 , and the improvement of two-photon over one-photon femtosecond photoassociation 
is expected to be less significant. Alkali and alkaline earth species where cooling and trapping 
of atoms is well established are considered. Carrier frequencies in the near-infrared are assumed 
throughout. 

A. Two-photon transitions in alkali dimers 

A one-photon transition cannot be controlled by interference between different photon pathways. 
A prerequisite for our scheme is therefore that no one-photon resonances be contained within the 
pulse spectral bandwidth. This excludes both potassium and rubidium as possible candidates. Of 
the remaining alkalis, only caesium shows a two-photon resonance between the ground state and a 
d-level in the near-infrared, cf. Table HI 

B. Two-photon transitions in alkaline earth dimers 

Of the alkaline earth atoms, calcium and strontium show two-photon resonances at the desired 
wavelength, cf. Table HI and in magnesium a near-IR four-photon transition to a d level is possible. 
In the following we will study calcium since potential curves for a large number of electronically 
excited states are known with sufficient accuracy^. ' The two-photon excitation scheme is shown 

* One may also consider transitions into potentials correlating to an S + P asymptote. While this is two-photon 
forbidden for atoms, it might be two-photon allowed at short internuclear distances. However, in that case 
photoassociation will occur only at distances below 20 Bohr. This corresponds to large detunings, and a two- 
photon control scheme is then not required. 

t We expect that our reasoning can be carried over to strontium with only slight modifications due to the similar 
electronic structure. 



in Fig. [2j Two atoms interacting via the X 1 S+-state are excited into bound levels of the (l) l H g - 
state or the (2) 1 E^"-state that both correlate to the 1 S + 1 D asymptote. Note that the 1 D term 
of calcium occurs energetically below the 1 P° term. The dipole coupling is provided by the states 
with one-photon allowed transitions from the ground state, i.e. the A 1 ^^, ^l' 1 ]!^, B 1 ^^ states 
and further electronic states at higher energies. The effective two-photon coupling, cf. Eq. ([5]), 
is determined by the dipole matrix elements of the one-photon allowed transitions and by the 
one-photon detunings. Compared to the atomic two-photon transitions studied in caesium^ and 
sodium 21 ! 22 ^ 25 , the one-photon detunings in the alkaline earths are much larger. Higher laser 
intensities will therefore be required. 

Fig. demonstrates the existence of long-range vibrational wavefunctions in a potential with 
1/i? 5 behavior for large internuclear distances R. Vibrational wavefunctions of the (l) 1 ]!^ excited 
state potential are plotted in Fig. EH for three different binding energies. The wavefunctions need 
to be compared to the initial scattering state, shown in Fig. [3b for two ground state atoms colliding 
with approximately 40 fiK: Large free-bound Franck-Condon factors are due to sufficient probability 
density at large R in the excited state vibrational wavefunctions. In the following population shall 
be transferred from the initial ground state atomic density into these vibrational wavefunctions 
close to the excited-state dissociation limit. The efficiency will be determined by the width of the 
dark atomic resonance, i.e. by how close to the dissociation limit excitation into molecular levels 
can occur. 

V. Two-photon photoassociation in calcium 

Photoassociation of two calcium atoms is studied by numerically solving the time-dependent 
Schrodinger equation for the effective two-photon Hamiltonian, Eq. Q, with a Chebychev 
propagator 29 . The Hamiltonian is represented on a grid using an adaptive grid step s i z o 30 i 31 i 32 . The 
two-photon couplings and dynamic Stark shifts are calculated from the atomic transition dipole 
matrix elements and frequencies found in Ref.—. The potentials are gathered from Ref. 28 . An 
initial state with a scattering energy corresponding to 40 ^K, cf. Fig. [3b, is chosen. 

Fig. H] presents calculations with transform-limited pulses of 100 fs full-width half-maximum 
and a central wavelength of 915 nm. They serve as reference point for the shaped pulse calcu- 
lations discussed below. Fig. EH shows the final excited state population including both atoms 
and molecules as a function of pulse energy. The arrows in Fig. 0H indicate pulse energies in the 
weak, intermediate and strong field regime that are employed in shaped pulse calculations. Rabi 



oscillations are observed as the pulse energy is increased. The contrast of the oscillations is reduced 
from 100% to 14%. This is a clear sign of the dynamic Stark shift for strong fields. For strong 
(intermediate) field, the maximum Stark shift, Eq. (|6j), amounts to 104 cm^ 1 (32cm -1 ) for the 
ground state and to 247 cm -1 (76 cm" 1 ) for the excited state. The weak field corresponds to a 
pulse energy of 0.94 /iJ and yields an excited state population of 5.0 x 10~ 3 , the intermediate field 
to 2.4 /iJ and P eX c{t final) = 2.9 x 10~ 2 , and the strong field to 7.6 and P eX c(t final) = 1-4 x 10 _1 . 

The projection of the final wave packet onto the molecular levels of the excited state, i.e. the 
excited state vibrational distribution after the pulse, is shown in Fig. 0b for transform-limited two- 
photon 7r and 2-7r-pulses. For a 7r-pulse, the population of the last bound level amounts to 1.8 x 10 -3 
compared to an atomic population of 1.4 x 10 _1 . For comparison, one-photon photoassociation 
with a lOps-pulse detuned by 4cm~ 1 achieves a population transfer into molecular levels on the 
order of 10 -4 which corresponds, depending on the trap conditions, to 1-10 molecules per pulseA 3 -^ 4 -. 
The arrows in Fig. 0b indicate the binding energies of a few weakly bound levels. Due to their 
smaller free-bound Franck-Condon factors, deeper molecular levels are populated significantly less 
than the last bound level. The minimum of the excited state vibrational distribution at v = —4 is 
due to the node of the scattering wavefunction at R ~ 68 ao, cf. Fig. [3b. 

Frequency-domain control is studied for weak and intermediate fields in Fig. The shaped 
pulses are obtained by applying a spectral phase function consisting of a 7r-step to a transform- 
limited lOOfs-pulse, cf. Eq. (|lip . The position of the 7r-step is varied. Atomic population transfer 
(black filled circles) can be strongly suppressed, by 10 -6 , for the weak field, see Fig. [5^. For 
the intermediate field, Fig. [5b, the dynamic Stark shift moves the levels during the pulse. The 
dark-pulse condition taking into account only the static resonance frequency, Eq. (jlOH . is then not 
sufficient anymore to strongly suppress atomic population transfer. The maximum suppression, 
i.e. the minimum of the final excited state population, is therefore reduced to be on the order of 
10 -3 . It is observed at the position of the 7r-step corresponding to Eq. (|12p . detuned by 42.2 cm -1 
from the carrier frequency for both weak and intermediate field. 

When population transfer into all molecular levels is considered, no difference between atoms 
and molecules is observed, compare the black filled circles and open green squares in Fig. [5] In 
this case the last two levels which have the largest free-bound Franck-Condon factors dominate 
the molecular population transfer. Since their binding energy is very small, less than 0.01 cm -1 , 
there is almost no detuning with respect to the atomic two-photon resonance, and the dark-state 
condition applies to the atomic transition and transitions into the last two bound levels alike. 
As the binding energy of the molecular levels and hence the detuning from the atomic resonance 



increases, the two-photon dark state condition which is defined at the atomic resonance becomes 
less applicable. Subsequently, molecular two-photon transitions become less dark, cf. open triangles 
and diamonds in Fig. [5j Atomic transitions are suppressed by about four orders of magnitude more 
then transitions into molecular levels with binding energies larger than 1 cm -1 for the weak field. 
This needs to be compared to the difference in atomic and molecular transition matrix elements. For 
weakly bound levels, it amounts to four to five orders of magnitude. The suppression of the atomic 
transition in Fig. implies therefore that about the same number of atoms and of molecules bound 
by more than 1 cm -1 will be excited. The excited atoms will be lost. However, this loss is sufficiently 
small, and the trap is not depleted. Due to the high repetition rate in femtosecond experiments, 
many photoassociation pulses can be applied before any significant depletion of the trap occurs. At 
intermediate field strength, the atomic transition relative to transitions into molecular levels with 
binding energies larger than 1 cm -1 is suppressed by only two orders of magnitude. In that case, 
loss of atoms from the trap becomes significant, and the control strategy of applying a spectral 
7r-step phase function is not sufficient to accumulate molecules by applying many photoassociation 
pulses. 

For intermediate and strong fields, pulse shaping in the time-domain is more appropriate since 
it allows to correct for the dynamic Stark shift. Time-domain control is studied in Fig. [6J The 
shaped pulses are obtained by chirping the pulse according to Eq. (fl6|) with the chirp time interval, 
At x , determined from Eq. (|17p . The final excited state population, normalized with respect to the 
transform limited case, is compared for atoms and molecules in Fig. [6] for intermediate (a) and 
strong (b) fields. It is plotted for increasing duration of the transform-limited pulse from which 
the chirped pulse is generated, keeping the pulse energy constant. For time-domain control, two- 
photon transitions for atoms are clearly more strongly suppressed than for molecules. The ratio of 
suppression of atoms versus molecules varies as a function of the transform-limited pulse duration 
between one and five orders of magnitude demonstrating a large amount of control. The overall 
increase of the difference between atoms and molecules with increasing pulse duration is rationalized 
in terms of the control strategy: The dynamic Stark shift together with the two-photon detuning 
leads to the accumulation of a phase in the atoms and molecules. It is countered by the laser phase 
<p(t), cf. Eq. Q. The difference between atoms and molecules becomes more perceptible for longer 
phase accumulation times. The modulations seem to be caused by a periodic accumulated phase of 
the instantaneous transition dipole. More analysis using local control theory is required. Overall, 
however, the suppression of the molecular transitions is much too strong to achieve any significant 
photoassociation yield. Analysing the molecular transitions for the weakly bound levels relevant 



in photoassociation, no dependence on binding energy is observed. Transitions into levels bound 
by 1 cm" 1 or more are equally suppressed as transitions into the last two bound levels. This is not 
surprising since the levels have a very similar dynamic Stark shift. Differences arise mostly from 
the different two-photon detunings, and they are much smaller than the dynamic Stark shift. It is 
hence not enough to enforce a minimization of the atomic transition probability. The molecular 
transition probabilities need to be maximized explicitly, e.g. employing local control and Eq. (I20p , 

VI. Summary and conclusions 

We have studied coherent control of femtosecond photoassociation employing two-photon tran- 
sitions. The goal of our study was to identify pulses which populate molecular levels close to the 
atomic transition without exciting the atoms themselves. An effective two-photon Hamiltonian 
was derived by extending previous work on atoms^ 3 - 1 ^ by the vibrational degree of freedom. In the 
weak-field to intermediate-field regime, frequency-domain control was pursued based on applying a 
spectral phase function. This leads to constructive or destructive interference between all pathways 
contributing to the two-photon absorptio n 19 ' 20 ' 21 ' 22 . Time-domain control is more appropriate for 
the intermediate-field to strong-field regime where the dynamics are dominated by dynamic Stark 
shift a 23 ' 24 . A temporal phase function allows to lock the atomic phase to that of the laser, keeping 
the atoms on or out of resonance. For both regimes, we derived analytical conditions to construct 
pulses which leave the atomic two-photon transition dark. These pulses were tested in a numerical 
study of photoassociation in calcium. An excited-state potential correlating to the 1 S+ 1 D asymp- 
tote was chosen. It shows a 1/fi 5 behavior at long range and supports weakly bound molecular 
levels with large free-bound Franck-Condon factors required for efficient photoassociation. 

In the weak-field regime, applying a 7r-phase step reduces the atomic transition probability by 
a factor of 10 -6 . Molecular levels with binding energies of about 1cm" 1 are sufficiently detuned 
from the atomic two-photon resonance that the dark condition does not apply. Transitions into 
these levels are suppressed by only a factor of 10" 2 . The resulting four orders of magnitude 
difference between atomic and molecular transitions is sufficient to counter the different excitation 
probabilities for atoms and molecular levels. It hence allows for employing femtosecond pulses 
without depleting the trap due to excitation of atoms. However, weak-field control implies that 
absolute excitation probabilities are very small. 

In the intermediate-field regime, the Stark shift becomes large enough to move the levels dur- 
ing the pulse. The spectral dark pulse condition taking into account only the static two-photon 



resonance is then not sufficient to suppress atomic transitions. A reduction of the atomic tran- 
sition probability by merely a factor of 10~ 3 was found. Applying a femtosecond pulse will then 
excite significantly more atoms than molecules. Since the excited atoms are lost, trap depletion is 
expected to set in after a few pulse cycles. 

Time-domain control based on maintaining a 27r-pulse condition for the atomic transition is 
applicable for both intermediate and strong fields. The ratio of suppression of atomic versus 
molecular transitions depends on the duration of the transform-limited pulses from which the 
shaped pulses are generated. It reaches up to five orders of magnitude. This finding is rationalized 
in terms of accumulation of different phases for atoms and molecules. In absolute numbers, however, 
the molecular transitions themselves are strongly suppressed. A significant photoassociation yield 
can therefore not be achieved by the simple control strategy based on the analytical result for 
the atomic two-photon excitation. Local control should be employed to maximize the molecular 
transition probabilities while keeping the atomic line dark. This is beyond the scope of the present 
study. 

In conclusion, an impressive amount of control yielding several orders of magnitude difference 
between atomic and molecular transitions could be demonstrated with both spectral and temporal 
phase shaping. However, the two-photon photoassociation efficiencies need to be improved in 
terms of absolute numbers of molecules per pulse to yield a feasible photoassociation scheme. 
In the strong-field regime, local control offers a means for explicit maximization of molecular 
transition probabilities. Moreover, our discussion could be extended from two-photon to three- 
photon transitions. Then, excited-state potentials correlating to S + P asymptotes and showing a 
1/R 3 behavior at long range could be employed. A 1/R 3 long-range potential supports spatially 
more extended vibrational wavefunctions than a 1/R 5 potential and exhibits an overall larger 
density of vibrational levels. Therefore, the efficiency of a shaped pulse which allows to excite 
weakly bound molecular levels in a three-photon transition while keeping the atomic resonance 
dark is expected to by far exceed that of the two-photon scheme. 
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atom 


transition 


wavelength 


Cs 


6s 7d 


2 x 768 nm 


Mg 


3s 2 -> 3s3d 


4 x 862 nm 


Ca 


4s 2 -» 3d4s 


2 x 915 nm 


Sr 


5s 2 -> 4d5s 


2 x 992 nm 



TABLE I: Two-photon transitions which are candidates for two-photon photoassociation in the near-infrared 
for alkali and alkaline earth atoms 



Fig. 1 In one-photon photoassociation, blocking the spectral amplitudes at and close to the 
atomic resonance prevents weakly bound molecular levels with large free-bound Franck- 
Condon factors to be excited (right-hand side). The majority of spectral components 
does not find any ground state population to be excited since the probability to find 
two colliding atoms at short internuclear distance is tiny (left-hand side, data shown 
here for Rb2). In femtosecond photoassociation employing a one-photon transition, the 
large majority of spectral components passes through the sample without any effect. 

Fig. 2 Two-photon photoassociation of two atoms colliding over the X l Yi£ ground state into 
the (l) 1 II 5 -state or the (2) 1 S+-state correlating to the s+d asymptote (solid lines). The 
dipole coupling is provided by the states with one-photon allowed transitions from the 
ground state such as the A 1 E^-state or the i? 1 S^-state (dashed lines). The example 
shown here is C&2- 

Fig. 3 Weakly bound vibrational wavefunctions of the (l) 1 !^ excited state for three different 
binding energies relevant for photoassociation (a) and scattering wavefunction of two 
ground state calcium atoms (b). 

Fig. 4 (a) The excited state population shows Rabi oscillations as the pulse energy is increased. 

The arrows indicate the pulsed energies which are employed for shaped pulses, (b) 
The projection of the final excited state wave packet onto the molecular levels, i.e. the 
excited state vibrational distribution after the pulse for a two-photon 7r-pulse and a 
two-photon 27r-pulse. The arrows indicate the range of binding energies relevant in the 
discussion of the results for shaped pulses. 



Fig. 5 Frequency-domain control with pulses shaped according to Eq. fjl 1[> : The excited state 
population of atoms (filled circles) and of molecules (open symbols), normalized with 
respect to the results obtained for a transform-limited pulse, is shown as a function of 
the 7r-step position for (a) weak and (b) intermediate fields. Population transfer can be 
strongly suppressed, by 10 -6 , for the weak field (a). For the intermediate field (b), the 
dynamic Stark shift sets in yielding a maximum suppression on the order of 10 -3 . When 
all molecular levels are considered, the last two levels (that have binding energies of less 
than 0.01cm -1 ) dominate and no difference between atoms and molecules is observed. 
For increasing binding energy of the molecular levels, and hence detuning from the 
atomic resonance, the two-photon transition becomes less dark. 



Fig. 6 Time-domain control with pulses shaped according to Eq. (j 16() : The excited state 
population of atoms (filled circles) and of molecules (open symbols), normalized with 
respect to the results obtained for a transform-limited pulse, is shown as a function 
of the duration of the corresponding transform-limited pulse for intermediate (a) and 
strong field (b). Two-photon transitions for atoms are more strongly suppressed than 
for molecules. However, the suppression of the molecular transitions is too strong to 
achieve any significant photoassociation yield. For weakly bound levels no difference 
in the suppression with increasing binding energy is observed, i.e. the levels have very 
similar dynamic Stark shifts. 
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FIG. 1: In one-photon photoassociation, blocking the spectral amplitudes at and close to the atomic reso- 
nance prevents weakly bound molecular levels with large free-bound Franck-Condon factors to be excited 
(right-hand side). The majority of spectral components does not find any ground state population to be 
excited since the probability to find two colliding atoms at short internuclear distance is tiny (left-hand side, 
data shown here for Rb2). In femtosecond photoassociation employing a one-photon transition, the large 
majority 
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FIG. 2: Two-photon photoassociation of two atoms colliding over the X l Y^ ground state into the (l) 1 !^- 
state or the (2) 1 E+-state correlating to the s + d asymptote (solid lines). The dipole coupling is provided 
by the states with one-photon allowed transitions from the ground state such as the A 1 S+-state or the 
B 1 S+-state (dashed lines). The example shown here is Ca2- 
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FIG. 3: Weakly bound vibrational wavefunctions of the (l) 1 ^ excited state for three different binding 
energies relevant for photoassociation (a) and scattering wavefunction of two ground state calcium atoms 
(b). 
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FIG. 4: (a) The excited state population shows Rabi oscillations as the pulse energy is increased. The arrows 
indicate the pulsed energies which are employed for shaped pulses, (b) The projection of the final excited 
state wave packet onto the molecular levels, i.e. the excited state vibrational distribution after the pulse for 
a two-photon 7r-pulse and a two-photon 27r-pulse. The arrows indicate the range of binding energies relevant 
in the discussion of the results for shaped pulses. 



1x10 
1x10" 



-3 1x10 



Oh 

H 



N 

13 



o 



1x10 
1x10 
1x10 
1x10 



1x10 



1x10 



1x10 



1x10 




' — • atoms 
°- -° molecules 

molecules 
with E u . ■ 

bind 

^>0.01 cm" 1 
l 



>0.1 cm 

>0.5 cm 

-i 

>1 cm 



l 



step position from resonance [ cm ] 



FIG. 5: Frequency-domain control with pulses shaped according to Eq. The excited state population 

of atoms (filled circles) and of molecules (open symbols), normalized with respect to the results obtained for 
a transform-limited pulse, is shown as a function of the 7r-step position for (a) weak and (b) intermediate 
fields. Population transfer can be strongly suppressed, by 10 -6 , for the weak field (a). For the intermediate 
field (b), the dynamic Stark shift sets in yielding a maximum suppression on the order of 10 -3 . When 
all molecular levels are considered, the last two levels (that have binding energies of less than 0.01 cm -1 ) 
dominate and no difference between atoms and molecules is observed. For increasing binding energy of the 
molecular levels, and hence detuning from the atomic resonance, the two-photon transition becomes less 
dark. 
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FIG. 6: Time-domain control with pulses shaped according to Eq. (fTH| : The excited state population of 
atoms (filled circles) and of molecules (open symbols) , normalized with respect to the results obtained for a 
transform-limited pulse, is shown as a function of the duration of the corresponding transform-limited pulse 
for intermediate (a) and strong field (b). Two-photon transitions for atoms are more strongly suppressed 
than for molecules. However, the suppression of the molecular transitions is too strong to achieve any 
significant photoassociation yield. For weakly bound levels no difference in the suppression with increasing 
binding energy is observed, i.e. the levels have very similar dynamic Stark shifts. 



